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We introduce bi-galileon theory, the generalisation of the single galileon model introduced by 
Nicolis et al. The theory contains two coupled scalar fields and is described by a Lagrangian that 
is invariant under Galilean shifts in those fields. This paper is the first of two, and focuses on 
the motivation and formulation of the theory. We show that the boundary effective theory of the 
cascading cosmology model corresponds to a bi-galileon theory in the decoupling limit, and argue 
that this is to be expected for co-dimension 2 braneworld models exhibiting infra-red modification 
of gravity. We then generalise this, by constructing the most general bi-galileon Lagrangian. By 
coupling one of the galileons to the energy-momentum tensor, we pitch this as a modified gravity 
theory in which the modifications to General Relativity are encoded in the dynamics of the two 
galileons. We initiate a study of phenomenology by looking at maximally symmetric vacua and 
their stability, developing elegant geometric techniques that trivially explain why some of the vacua 
C ' have to be unstable in certain cases (eg DGP). A detailed study of phenomenology appears in our 

companion paper. 
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I. INTRODUCTION 



Oj' When Urbain Le Verrier first noticed that the slow precession of Mercury's orbit could not be explained 

QJ I using perturbative Newtonian gravity and the known planets, he postulated the existence of a hitherto 

1-^5 ■ unseen dark planet, which he called Vulcan [l|. Of course, thanks to Einstein, we now know that what was 

really required was a new theory of gravity, specifically, the General Theory of Relativity [2J. 150 years later, 

tH- . flat rotations curves of galaxies [3] and accelerated cosmic expansion Q have led most of us to follow Le 

J> ' Verrier's approach and to postulate the existence of dark matter [5^ and dark energy [6] respectively. For 

^^ , dark matter, bullet cluster observations suggest that Le Verrier's approach is probably the right one \7]. 

Cn ■ The Le Verrier approach to dark energy is less compelling. Whereas extensions of the Standard Model offer 

some very good dark matter candidates, they have been nothing short of hopeless in offering a satisfactory 
candidate for dark energy. It is certainly worth taking a lesson from history and exploring the possibility 

C^^ ' that a new theory of gravity may be required to explain the observed cosmic acceleration. 

^^ , In order to play a role in the accelerated expansion, we need gravity to deviate from General Relativity 

^-^ ' (GR) at large scales, and in particular at the scale of the cosmic horizon Hq^ ~ lO^^m. Such deviations 

should not be felt all the way down to arbitrarily short distances because General Relativity is very successful 
in describing classical gravity up to the scale of the Solar System, at about lO^'^^m. Modifications of General 

k><( , Relativity are most easily understood at the level of linearised theory. In GR, the gravitational force is 

j_] ■ mediated by a massless spin 2 particle carrying two polarisation degrees of freedom. In modified gravity we 

C^ ' expect additional degrees of freedom to be present. For example, in Pauli-Fierz massive gravity the graviton 

carries 5, rather than 2, polarisation degrees of freedom [8[. Clearly if they are to play a role in accelerated 
expansion without spoiling the phenomenological successes of GR at shorter distances, the additional degrees 
of freedom should be suppressed up to Solar System scales, but not on cosmological scales. This should 
occur without introducing any new pathological modes. In massive gravity, the longitudinal graviton mode 
is problematic since it gives rise to the so called vDVZ discontinuity in the graviton propagator, even at 
scales where the graviton mass is negligible Q. However, this mode can be "screened" on Solar System scales 
thanks to non-linear interactions and the Vainshtein mechanism [lO| , but only at the expense of introducing 
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a new pathology in the form of the Boulware-Deser ghost [Tl|, 113 • 

Extra dimensions and the braneworld paradigm |13| hav e op ened up some intrig uing possibihties for 
modifying gravity at large distances (see, for example, |14| - [l7| ). The DGP model [ij] is probably the 
most celebrated example of this. The model admits two distinct sectors: the normal branch and the self- 
accelerating branch. The latter generated plenty of interest since it gave rise to cosmic acceleration without 
the need for dark energy [1^. However, fluctuations about the self-accelerating vacuum were found to suffer 
from ghost-like instabilities, ruling out this sector of the theory [1^ [20|. Although the normal branch is 
less interesting phenomenologically, it is fundamentally more healthy and is the closest thing we have to a 
consistent non- linear completion of massive gravity. Here the graviton is a resonance of finite width Hq^ as 
opposed to a massive field. At short distances, r <^ H^^ the brane dynamics does not feel the width of 
the resonance and the theory resembles 41? GR. At large distances, r > Hq^, however, the theory becomes 
five dimensional as the resonance effectively decays into continuum Kaluza-Klein modes. The Vainshtein 
mechanism works well on the normal branch of DGP, screeniiig the longitudinal graviton without introducing 
any new pathological modes, in contrast to massive gravity |12l |21|. 

Much of the interesting dynamics is contained in the properties of the longitudinal scalar, tt. One can 
identify this mode with the brane bending mode. As this mode only excites the extrinsic curvature of the 
brane, and not the brane or bulk geometry, it becomes strongly coupled well below the 4Z) Planck scale |20| . 
It is precisely this strong coupling that allows the Vainshtein mechanism to take effect as non-linearities 
become iinportant at larger than expected distances, helping to shut down the scalar on Solar System 
scales [l^, [2J] . We can isolate the dynamics of the scalar by taking the so-called decoupling limit around flat 
space, in which the iD and 5D Planck scales go to infinity, but the strong coupling scale is held fixed. In 
this limit the tensor mode, h^^^, decouples from the scalar, tt, to all orders |20l. |22||. The 7r-Lagrangian now 
contains a plethora of information about both branches of solution, which can be studied much more easily 
than in the full theory. 

Now an important feature of the 7r-Lagrangian is that it exhibits galileon invariance [24[. That is, the 
action is invariant under tt — >■ tt -I- 6^x^ + c, where 6^ and c are constant. This essentially follows from the 
fact that the brane bending contribution to the extrinsic curvature is given by Kfj^^, ~ dfj^d^-K. Equivalcntly, 
Poincarc transformations in the bulk should not affect the physics of the brane bending mode. If the brane 
position is given by y = 7r(x), then after an infinitesimal Poincare transformation, x" — > a;° -I- e°'bX^ -f v"' 
the brane position is given by y ~ 7r(x) -I- e^^cc'' -|- v^ . Since the physics of the brane bending mode should 
remain unaffected it is clear that it should exhibit galileon invariance. 

Motivated by the DGP model, Nicolis et al [2J1 developed the idea of galileon invariance and wrote down 
the most general Lagrangian for a single galileon degree of freedom (see [25| for some early work along similar 
lines). Remarkably the theory contains only five free parameters, enabling the authors, and others J26|, to 
complete a thorough analysis. Although the Galilean symmetry is broken when we include coupling to 



gravity [27[ , covariant completions of the general galileon action can be identified with the position modulus 
of a probe DBI brane 28]. In any event, it seems clear that any co-dimension one brane model with large 
distance deviations from GR, will indeed admit a galileon description around flat space, in some appropriate 
limit. This follows from the discussion of the brane bending mode in the previous paragraph, and from [23j . 
where it is argued that any viable infra-red modification of gravity leads to strong coupling, and is therefore 
expected to admit a non-trivial decoupling limit. 

Recently there has been plenty of interest in branes of higher co-dimension and the role they might play 
in modifying gravity at cosmological scales |29l- [38| . An example of this is the cascading DGP model [295, ^^ 
which one has DGP branes within DGP branes |29l- [3l |. In the 61? version one has a co-dimension 2 DGP 
3-brane intersecting a co-dimension one DGP 4-brane [29| . Higher co-dimension brane models seem to have 
extra features, such as degravitation [30,[33|, or perhaps even self-tuning |29l [33 - l34| . in which the geometry 
of the 3-brane is insensitive to a large 3-brane vacuum energy. Now to specify the position of a co-dimension- 
N brane we need to specify the values of N "bulk" components: yi — 7ri(a:), j/2 = t^2{x), . . . ,yN — 'rr]sf{x). 
For a viable model giving rise to large distance deviations from GR, there should be strong coupling and 
therefore an appropriate limit in which the brane bending modes decouple from the tensor mode around flat 
space. Then the physics of each brane bending mode should be unaffected by Poincare transformations in 
the bulk, and we conclude that we are left with an A^-galileon theory. 

In this paper we will focus on the case of two coupled galileons, which we will call bi-galileon theory. Of 
course, this has particular relevance to co-dimension 2 brane models. Indeed, in section [TTl we will show 



explicitly that the 5D cascading cosmology model is one example of a bi-galileon theory. This model is 
derived from the original 61? cascading gravity model [29|, |30[ by computing the boundary effective field 
theory in 5 dimensions, and taking a suitable decoupling limit. Although a detailed study of phenomenology 
will appear elsewhere [4l|, this paper is devoted to the motivation and formulation of the most general 
bi-galileon theory. In section IIIII we will pitch our theory as a modification of gravity by coupling one of 
the galileon fields to the energy-momentum tensor. We will not couple the graviton directly to the scalars, 
which amounts to ignoring their backreaction on to the geometry, as was done for the case of the single 
galileon [2J|. In section HVl we will construct the general bi-galileon theory. We will also state the extension 
to any number of coupled galileon fields. We initiate a study of phenomenology in section |V] by looking at 
maximally symmetric vacua and their stability. Section [VII contains a half-time analysis. 



II. EXAMPLE: CASCADING COSMOLOGY 



There have been a number of attempts to extend the DGP model to higher dimensions. Although the 
earliest of these ran into problems with ghosts and instabilities [401 j the recently developed cascading mod- 
els [29I - I3H have been more successful. These models exhibit degravitation and/or self-tuning, and seem to 
be free from ghosts, at least if the 3-brane tension is large enough. In this section we will compute the AD 
boundary effective theory for the cascading cosmology model [31j] . By taking the appropriate limit we will 
decouple the scalar sector and show that it corresponds to a bi-galileon theory. 

The cascading cosmology model [31[ is closely related to the 6D cascading DGP model ^^. The latter 
contains a DGP 3-brane sitting on a DGP 4-brane in 6 dimensions. The action is given by 
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If we have a hierarchy of scales Mq < M^ < A/4, then the Newtonian potential cascades from 1/r to 1/r^ to 
1/r^ as we move further away from the source. The behaviour crosses over from 4D to 51? at a scale m'^^ 



and from 51? to 6D at a scale 



where 
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In the original model [29] both crossover scales are taken to be of order the cosmic horizon, so that ms ~ 
TOg ~ i/o ~ 10"^^ eV. Taking M4 ~ Mpi - lO^^ GeV, it follows that M5 ~ 10 MeV and Mg - meV. Note 
that bulk quantum gravity corrections are expected to kick in at Mq^ ^0.1 mm. 

Now, one can integrate out the bulk, just as in [20], to obtain the effective 51? theory. It turns out that 
the longitudinal mode of the 51? graviton becomes strongly coupled at a scale 
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or, equivalently, at distances of order Ag ~ 10^ km. Now at distances much shorter than the cosmic horizon, 
but much larger than the millimeter scale, we can happily take MsjAfg — > 00, whilst holding Ag fixed. This 
corresponds to the decoupling limit in the effective 5D theory [29|, |3i|. The cascading cosmology model [3i| 
represents a non- linear completion of this decoupled effective theory^ . The action is given by 
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^ This non-linear completion is certainly not unique but one can reasonably expect it capture the salient features of the original 
6D theory. 



where the fuh 51? spacetime is made up of two copies of A4, identified across a common boundary, dAi, 
given by the 3-brane. Here gabdx°'dx^ and q^^dx^dx'^ are the bulk and brane metric respectively, with 
corresponding Ricci scalars R^ and R^. The extrinsic curvature of the brane is given by K^i, = \Cnq^_Lv 
where £„ is the Lie derivative with respect to the inward pointing normal, n"' , to dM. The bulk scalar ■k{x) 
is the bD longitudinal mode and the parameter a — -j^. 

If we choose coordinates x"" = (y, x''), so that the brane is at y = and the bulk extends into y > 0, then 
it is convenient to write the bulk metric in ADM coordinates 

ds^ = gabdx'^dx'' = N^dy^ + <j^^(dT^ + N''dy){dx'' + N^dy) , (6) 

where N is the lapse function and N'^ is the shift vector. Then 

K^. = 2^ {dyq^. - 2D^^N,^) , -C^tt = 1 {dy - N^D,) it , (7) 

where D^ is the covariant derivative with respect to the brane metric g^i/. 

To calculate the boundary effective action, we essentially follow steps completely analogous to the ones 
taken in [20] ■ In order to integrate out the bulk fields we must choose a gauge. We therefore set 

gab = Vab + hab , <7^zy = VfJ.1^ + hf_,u , (8) 

and add a bulk gauge fixing term 

Af3 /• s S 

Sbuik,yf = -^ /^ ^'^ (P- + 2^a7r){F'^ + ^d'^n) , (9) 

where 

F- = dbh-' - ^d'^h , h = Tj'^'hab ■ (10) 

Classically this imposes the gauge Fa — — |(9a7r. To quadratic order we have 
^2^ = Ml fd^xdy i ("/i"^ - ^hrjA D5hab + ^nO^h + ^TrDg^ 

+ Ml f d^X i {h^''' - ^hil^'"') dyh^, + ^^N'^dyN^ + -^^^dyh^ - ^{hyydyh^ + h^dyhyy) 
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+ -^hyydyhyy + -TTdykyy + -TldyTT + N^ ( /^ ^ ^O'^kyy + -8'' IT 

+ ^ /_o '^''^ i (^""^ " ^''"^''0 °'^^'' ^ ^^''^^ ' ^^^^ 



r = d,h^^- - -d^K , hi = V'^'^h^, , (12) 



where 

1 

2' 
and we have re-defined iV^ = h^y . The bulk equations of motion 

Ds/lab-O, DsTT^O (13) 

have solutions 

hab^ e^^^hab\y=0, TT = e^^'^7r|y=o , (14) 



where A = -y/— 04. Now we integrate out the bulk and fix the residual gauge freedom left on the brane by 
adding a brane gauge fixing term, 

Sbrane,gf = " "^ / d^X (/^ + m^N^X^ + m^N^^) . (15) 

4 Jm 



This gauge fixing eliminates the mixing between N^ and /^ and we are left with 
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We can diagonalise this action by making the following field re-definitions 



(16) 



hf^i, = hf^i, + (tt + 0)?7pi. , N^, = N^,^ d^{TT + 



TO5 



hyy = i'K (A + m5)(7r + 0). (17) 



This gives 
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\hiV''"j (04 - TO5A) V - ^™5iV^(A + m5)7V^ 
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-(/)(n4 - m5A)(/i ■7r(n4 + 2m5A)7r . 



(18) 



Note that the tt field is a ghost. This can be avoided by increasing the vacuum energy on the brane [29|. In 
addition, to quadratic order, the matter coupling is given by 
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52S.maUer = X / d^'xh^^T^^" = - / d'^X h^^T^" + (tT + (/.)T . 

^ Jy=Vi ^ •Jy=0 



(19) 



We now consider higher order terms. Schematically it is clear from Eq. [5]that these will take the form 

d^x Af|A(7r)«i (/i^,)^i {N^r^ {hyyf' , (20) 

d^x MlaA^inr- ( V)'^ (^p)'^ (/ij.j.)''^ , (21) 

d-'x Mjls'ih^.p , (22) 



where Ui, . . . ,di are non- negative integers satisfying ai + bi + Ci + di > 3. Note that we must have 02 > 1 
since terms of the form Eq. [5T]stem from the boundary n term given in Eq. [5l After canonically normalising 
the fields, it is easy enough to check that the largest contribution comes from a term like Eq. [20l with 
ai = bi = 0, ci + di = 3. This term is exactly like the cubic term computed in [20|, and comes from the 
following term in the action ([5]) 



Mi 
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Since we want ai = 0, we can switch off e"'^'^'^, and write this term in ADM form. 

Ml I d'^xdy^f^N{Ri + K'^ - K^.^K^") . 



(23) 



(24) 



Of course, we also want 61 = 0, so we need to switch off h^^. This means q^u — rj^^y, and so i?4 = and 
K^y — ~jrd(^p^Ni,y Given that A^ sa 1 + hhyy, the cubic order interaction term is given by 



-^ I d^Xdyhyy [{d^N^^f - d^^Ny^d'^N'^ 



(25) 



In the limit where m^ — >■ 0, the strongly interacting mode can be parametrized in the bulk as 

A^^-a^*, hyy^2dy'i, TT ^ , (26) 



/ipy ^ , 



where 05^ = 0. This satisfies the bulk equations of motion, as well as the gauge choice Fa = — |9a7r. Indeed, 
the strongly interacting mode corresponds to a trivial bulk geometry, as expected [2Q] . Furthermore, we see 
from Eq. [T7]that the boundary values for the fields go like 7V^ ~ ^9p(7r + 0), and hyy ~ — ^A(7r + 0) so 

we take ^ — e^^^ I ^^i^ J . Eq. [25]can now be written as 

-Ml I d^xdy dy^ [(04*)' - {d,,d,^f] = ^ fd^xdy dy [d^.^d^'^ai^ 

\ d^xd^,{n + <t>)d^{Ti + (t>)Ui{Ti + 0) . (27) 



2m5_ 

Now if we canonically normalise, so that tt = a/§ -p- and (f) — \/§-^, then we see that the cubic interaction 
term goes like 

^ j d^xd^,{n + i)d^'{T^ + ^)ni{TT + ^), (28) 



where 

As = {mlMif/^ - 10""cV . (29) 

This represents the scale at which the AD scalars becomes strongly coupled, and corresponds to a distance 
Ag"^ ^ 1000km. We can now take the formal limit^ 

M4,M5,Me,Tfj_,y ^ 00 , A5 = const , Ag = const, — ^ = const . (30) 

Note that this requires that TO5, mg — ?> 0. It is easy to see that the leading order matter coupling remains. In 
contrast, the large interaction term Eq. [25]is the only interaction term to survive this limit. As a result, the 
tensor h^^i, and vector Nf^ decouple from one another and from the scalars n and (j) to all orders. In terms of 
the canonically normalised fields, the action is given by 

^) ^M^ + ^)d^(n + ^)n,in + ^) + V^ + ;^(^ + ^)|^ ' (31) 

where h^^, = ^h^^^ and iSf^ — M^ N^^. Focussing on the scalar sector, the relevant piece of the action is 
given by 

S..^ = j d'x ^4>n,$ - ^nO^n - (^^) -1 5^(^ + 0)9^(71 + ^)n4(7r + 0) + -L(^ + <^)-^ . (32) 

The resulting equations of motion are explicitly Galilean invariant, 
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^ In taking this limit we inevitably project out the pure DGP scenario, since that would require us to take Mg = 0,Ag = 0, 
which is, of course incompatible with Eq. 1301 



That is, they are invariant under 

^ ^ TT + 6^.T^ + C , (j)^ (t> + &^.X^ + c , (35) 

for constant 6p,6^,c and c [i24i |. Thus we confirm what we expected. Cascading cosmology admits a bi- 
gahleon description in the full decoupling limit. The limit corresponds to taking the Planck masses to 
infinity but holding the 4Z5 and 5D strong coupling scales fixed. The latter correspond to 1000 km and 
1000000 km respectively, whilst the Planck lengths in the original model are all submillimeter. In taking 
the decoupling limit we have sent all the next largest interactions to zero, so our description is only valid at 
scales where these additional interactions may be ignored, up to the horizon scale. These interactions may 
have an important role to play in understanding various aspect of the full theory (eg. Vainshtein effects) , 
but since our goal was merely to demonstrate that a bi-galileon description exists in principle at certain 
scales, we do not concern ourselves with that here. 

The TT, 0-Lagrangian associated with Eq.[32]is the analogue of the 7r-Lagrangian in DGP gravity [20,[23. It 
should capture local physics on the 3-brane, corresponding to normalisable fluctuations in the 51? bulk. Note 
that this means it cannot capture the self tuning behaviour of the full 6D theory. This is because the self 
tuning solutions correspond to non-local, non-normalisable fiuctuations about the vacuum of the effective 
5D theory described by the cascading cosmology model ([S|). We could, of course, consider normalisable 
fluctuations about a non-trivial self tuning solution. We expect that this would also give a bi-galileon theory 
in the appropriate decouplin g lim it. For sufficiently large 3-brane tension it would actually correspond to a 
ghost-free bi-galileon theory [2^, unlike the example given here. However, a detailed analysis is outside of 
the scope of this paper. 

III. BI-GALILEON MODIFICATION OF GRAVITY 

The cascading cosmology model described in the previous section is an example of a modification of 
General Relativity by two additional galileon fields, tt and (p. In a generic theory of modified gravity, the 
amplitude for the exchange of one graviton between two conserved sources, T^i/ and T', is given by 

A = T^.D^^-'^Pt'^p , (36) 

where D'^'^"^ is the graviton propagator. Now, for General Relativity, the amplitude is given by 

where T = T^^. We will be interested in the case where gravity gets modified by two additional scalars so 
that locally we have 

SA^A- Agr. = T^—r + T^—T' . (38) 

q;iU4 a2LJ4 

Such a theory can be described by the following action 

d'^x —^yf^R{g) ~ — 7rin47ri - — 7r2n47r2 + -h^uT'"' + niT + TTaT -|- interactions , (39) 

where 5^1, — rjf^i, + h^j^^. The fiuctuation /i^^ is identified with the graviton in General Relativity since it 
satisfies the same linearised field equations. Therefore, for a given source and boundary conditions, the 
solution h^u coincides with the linearised GR solution. This statement is true to all orders in the limit, 

Mpi, ttj, T^"" ^- 00 , —4- = const , — — = const , (40) 

M , Mpi 



where the scalars decouple from h^i,. 
fluctuation 



The physical graviton differs from its GR counterpart by a scalar 



hf_,„ = /ipy + 2(7ri +7r2)77, 



pi' > 



(41) 



since the physical metric seen by matter is given by g^^, — rj^^i, + h^^ . 

Let us consider the decoupling limit pO)) . with the additional assumption that the strength of all scalar 
interactions are held fixed. This basically means we are neglecting the backreaction of the scalars on the 
geometry, so we can consider them as fields on Minkowski space. Note that we retain scalar self-interactions 
because they are essential for bi-galileon corrections to GR to be extrapolated from O(10~^) in the solar sys- 
tem to 0(1) at cosmological scales. Significant modifications to gravity at Hubble distances would inevitably 
have a role to play in the dark energy problem. The action is given by 
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hhrj^i, ) -|- ... is the linearised Einstein tensor, and the Lagrangian Utt^.-k^ only 



depends on the scalars tti and t^i. As we can see from Ea. l31[ the action for the cascading cosmology certainly 
takes this form in the decoupling limit. In Eq. [3I]the scalar action is Galilean invariant. To generalise this, 
we are interested in bi-galileon modifications of gravity for which £711, 7r2 is the Lagrangian for a bi-galileon 



theory. This means that if tt^ — ^ tt^ -I- ihi)^ 
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total derivative. 



Given the form of the physical graviton (|4ip , it is convenient to define new galileon fields tt = tti 
£, = TTi — TT2, SO that only one of the scalars couples to matter. The action now takes the form 
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where C-j^,^ corresponds to a bi-galileon theory. The physical graviton is now 



7r2 and 
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Although only n couples to matter directly, S, couples to matter mdirectly, through its mixing with tt. 
possible for the two scalars to be mixed at quadratic order, as well as at higher order. 

Given that we have taken all our fields to propagate on Minkowski space, one might be forgiven for thinking 
that we cannot say anything about cosmology, or perturbations about cosmological solutions, including de 
Sitter vacua. Fortunately, this is not the case. We can think of any spatial flat cosmological spacetime as a 
local perturbation around Minkowski space 2j|. In this context, local means local in both space and time, 
at sub- horizon distances and over a sub- Hubble time. If we take here to be x = and now to be i = 0, then 
for If I < iJ-i and \t\ < H'^, we have [H 



ds^ 



-dr^ + airfdf 
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-H^\x\^ + -{2H + H^)t'^ 



{-dt^ + dx^) 
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where the Hubble scale H and its time derivative H may be evaluated at i = 0. This is a perturbation about 
Minkowski space written in Newtonian gauge 
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where the Newtonian potentials are 

$ = - 



-H^\x\'^ + -{2H + H^)t^ , * = -$. 
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For a given cosmological fluid, the corresponding GR solutions have Hubble parameter Han. Since h^i, 
agrees with the linearised GR solution, we have hu = —2^cR, hij — "^^GRSij, where 



*Gi? = --H, 
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Now in our modified theory the physical Hubble parameter (associated with /i^,y) differs from the corre- 
sponding GR value, H ^ Hgr- Due to Eq. |44l we have a non-trivial scalar 

TT = $ - ^GR ■ (49) 

The other scalar, ^, follows from the equations of motion. Note that a Galilean transformation vr — ;> tt -|- 
b^x^ + c merely corresponds to a coordinate transformation x^ ^ x^ — cx^ + ^{x^^x'^b^ — 2b,jX^x^) in the 
physical metric. 

Of particular interest are maximally symmetric vacua, for which H and i^cfl. ^-re constant. Indeed, if the 
vacuum energy is a, we have Hqj^ — cr/3Mpj. It is easy to check that one scalar now takes the particularly 
simple form tt — —jk^^Xf^x^ where fc^r — H^ ~ ^gr- ^^ principle the other scalar can be anything, but we 
will typically assume that it takes the same form on maximally symmetric solutions, so that ^ — —^k^x^x^. 

Ultimately, we would like to promote our galileon theory to a fully covariant theory of gravity, coupled 
to a pair of scalar fields. Whilst this will inevitably break the Galilean invariance [23|, we should expect 
to recover the galileon description in the decoupling limit Mpi — > oo. One does have to be slightly careful 
not to introduce any higher derivatives by accident when covariantising the theory, but this can be done 
using additional non- minimal couplings [27]. As a step towards this, let us first transform the action (j43p to 
"Jordan" frame by rewriting everything in terms of the physical metric (j44p . This gives 

d^x - -^hf^^Sh^, - M^iTrid^.d.h^" - d^h) + -h^.T^" + £,,5, (50) 

where ^C^r,^ = C„.^ — iAAnd^TT. The covariant completion of this theory will take the form 



^cov — / " X\/ g 



— — (1 — 2ti)R -t- £mattcr + -^scalar 



(51) 



where >Cmattcr[5; ^n] is the matter Lagrangian, describing matter fields, ^„, minimally coupled to gravity. 
The scalar contribution, £scaiar[5; 7r,C]: is constructed as follows: take £,1.^ and let rj^i, — >■ 5^11,, d^ — J- V^, 
then add non-minimal couplings to curvature to compensate for the higher derivative terms in the equations 
of motion. 

Actually, this completion is useful in that it enables us to establish the validity of our galileon description. 
Recall that this neglects the backreaction of the scalars onto the geometry. From Eq. [nH we see that the 
scalars contribute an energy-momentum tensor 

^s^calarlS; 7^, C] = -J= J— / d'^X^/^ £,ealar (52) 

V y "HtJ-v J 

The gravitational field equations arising from the covariant action (|5ip are given by 

M^ 1 

^(1 - 2^)0^"^ -I- M2 (V^V-^ - g^'^U)!, - {T^^ + T^^J - (53) 



Although the galileon description includes non-linear effects at the level of the scalar field equations, it 
neglects them at the level of the gravity equation. Expanding Eq. [53]about Minkowski space, g^,y ~ rifj_^+hf^^, 
we find that 

M"^ 1 1 

-^Eh'^'^ + M^ (d'^d'^ - rj^'^d^) TT - -r^'^[r;; *„] = -TZuAv, ^, C] + • • • (54) 

where ". . ." denotes terms that are 0{h) suppressed relative to the terms on the LHS of Eq. [Ml The 
gravity equation in the galileon description corresponds to neglecting everything on the RHS of Eq. [Ml 
When is this justified? For a small graviton fiuctuation, we can consistently neglect those the terms that are 
0{h) suppressed, denoted by ". . .". However, even then, the backreaction of the scalars on to the geometry 
can be significant. In other words, it is not clear that the energy momentum of the scalars, evaluated on 
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Minkowski, can be consistently neglected. To justify neglecting this backreaction we conservatively require 
thatr,';:,,Jry,;7r,^]«M2£:/,M^ 

Provided this condition holds, we will see that our bi-galileon theory allows for a rich and interesting 
phenomenology, to be studied in detail in our companion paper 41]. As we have just discussed, we can 
describe cosmological solutions, such as de Sitter space, as local perturbations about Minkowski. We will also 
be interested in fluctuations about these solutions, and in particular about non-trivial maximally symmetric 
vacua. These just correspond to fluctuations in the galileon fields about their background values. As we 
mentioned earlier, the non-linearity of the galileon action will be essential to the success of escaping the local 
gravity tests and producing interesting cosmological effects simultaneously. 

IV. CONSTRUCTING THE BI-GALILEON THEORY 

Let us now construct the most general bi-galileon theory. Our task is to find the most general Lagrangian 
describing two scalars, /Itt,^, that is invariant under the Galilean transformation 



TT + bnX^ + C , 



£.^S.+ bf.x'' + c , 



(55) 



This amounts to requiring that the corresponding equations of motion are built exclusively out of second 
derivative terms d^d^TT, dadjs^ contracted with the background Minkowski metric . Single, or zero derivative 
contributions to the equations of motion break the Galilean invariance, whereas higher derivatives introduce 
extra degrees of freedom and ghosts. 

The Galilean symmetry greatly constrains the structure of the Lagrangian. For example, under an in- 
finitesimal constant shift tt — >■ tt -f c, the Lagrangian C^r.^ varies as 



-^TT.^ ^ ^7T,^ 






(56) 



However, as this shift happens to be a Galilean transformation, we know that any change in the Lagrangian 



must be a total derivative. It follows that dCn^^^/dTT 
goes like 



duJ'^, for some J^ . Now the tt equation of motion 



— / d'^x Ct^.i 









(57) 



which is clearly a total derivative d^ [J''' — 9£^^^/9((9p7r) + ...]. Indeed, the n equation of motion is a total 
derivative at each order in tt and ^. We can apply a similar argument to the i^ equation of motion. 

We can now formulate the problem as follows: at {m, n)-th order (rri-th order in n and n-th order in ^), 
what is the most general total derivative built out of second derivative terms? One obvious total derivative 
is given by 



£rr.,n^{m + n)%\...S^^ 



§P1 



■■CAd.^d-'^TT) . . . {d,^d''-n) (ap,a-o . • . (5p„a-o • 



(58) 



It is easy to see all such terms with m + n > A vanish in four dimensions. We will now write out 
all the non-vanishing £m,n explicitly in 4D. For notational convenience, we suppress the Lorentz indices 
and separate Lorentz scalars by '•' when it is confusing. For example, we have ddirdd^, ■ ddirdd^ = 
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{dfj,d,^n){d'^ d^£,){dadj3Tr){d^ d°'S,) ■ For m> n, the non-vanishing terms are 
Si^o = On, 

£-2,0 = (n7r)2 - (aa^)2 , 

£i4 = nnn^ ~ ddnddC , 

£3,0 = (n7r)3 - 3n7r(9a7r)2 + 2(aa7r)3 , 

£■2,1 = (n7r)2nC - n^(o»a7r)2 - 2n7r957ra9C + 2(997r)2aaC , 

£:4,o = (^7r)''-6(^7r)2(aa7^)2+8^7^(aa7^)3 + 3((9a7^)2•(aa7r)2~6(a97^)^ 

^3,1 = (□7r)3n^ - iOTrn^iddTT)^ ~ 3{aTT)^ddTTdd^ + 2a^{ddTTf 

+6aTT{ddTTfdd^ + 3(997r)2 • ddnddS, - 6{ddTrfdd^ , 
£2^2 = (□7r)2(n02 - (nO^(a57r)2 - (Dnfidd^)^ - AnTiD^ddndd^ 

+4a^{ddTrfdd^ + ADTTddTTiddO'^ + {ddnf ■ {dd^f 

+2ddTTddi ■ ddnddi - A{ddnf{dd£,f - 2{ddndd£,f , (59) 

For m < n, we take £m,n — '^n.mU^^?- Note that in each case £rn,n — (□7r)'"(n^)" + 

We win now prove that Eq. [SBlis the only total derivative of order (m, n) that is second order in spacetime 
derivatives. The proof works along the same hnes as the corresponding single galileon one given in 24]. 
First, assume there is another suitable total derivative £,'„„. Now consider £,'„ „ as a Lagrangian density. 
As £,'„ „ is a total derivative, the resulting vr and ^ equations of motion must be trivial 

|- J d'x C„ = , ^ / ^'^ C,„ = . (60) 

Without loss of generality, a general term in £^^ „ goes like^ 

(DTr)'" (nO^° (ddny^ {dd^^' • . . . • {ddTry" {dd^^" , (61) 

where io + . . . -\- ik = m, jo -\- . . . + jk = n- Let us now consider its contribution to the resulting equations of 
motion (j60p . A term in the tt equations of motion can be obtained by shifting dd (along with the Lorentz 
indices) from one ddir to another ddn or another ddS, and then dropping the n without derivatives. Therefore 
Eq. ini] must contribute a term proportional to 

{Dny^inO'^'DddTriddny^iddO'' ■■■■■ {ddny" {ddO^" , (62) 

in the equation of motion. As the overall equation of motion is trivial, this term must be cancelled by the 
variation of another term in £^ „. Aside from Ea. l61[ there is only one other possible term that can produce 
a term like Eq. IHH in the tt equation of motion. It is given by 

{D7ry'>+\no^° {ddny'-\ddO^' ■■■■■ {ddTry- {ddiy- . (63) 

It follows that this must also be in 5^ „. One can easily imagine a similar argument applies to ^. Now, if 
we apply these arguments iteratively, we would conclude that a term proportional to (n7r)™(n^)" must be 
in £1^ „. Then we immediately run into a contradiction, because linearly combining £m,n and f^ „ should 
produce another total derivative which is free of (n7r)™(n^)". This proves that the only suitable total 
derivative at order {m,n) is £,„,«• 



^ Note that we have suppressed permutations of tt and ^ within Lorentz scalars; for instance, we collectively list the Lorentz 
scalars dd-Kddidd-Kddi and {997r)2 (99^)2 as {ddTrfiddS)'^ 
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Returning to the full theory, the equations of motion for n and ^ in general correspond to linear combina- 
tions of the £m.n'- 

— tax /-^TT.f — / ^m.n^m.n ; 

On I ^-^ 

— / d*X £„^^ = ^ bm^nSm,n ■ (64) 

Since these equations are derived from the same Lagrangian, Ct^^^, we anticipate an integrability condition 
relating the coefficients am,n and bm,n- To construct the Lagrangian, we note that equations of motion with 
the correct form can be derived from the following 



^TT,^ — / ^ L'm.n J (,"&] 



m+n^A 



where 



To see this note that 



(ttm^nTT + Pm,n£,)^m,n ■ (66) 



6_ 

which yields the following 



— I d X i^ni,n ^^ \^^ I i- )(^m.7i^m.7i \ ?^Pm,nt^m— l,n+l 



A r 

— d'^xC^^^^ ^ (TO+l)(am,«+/?m+l,«-l)'?m,n (A,-l = 0) , (67) 

•^ n ,-^ ^, I — -^ /I 

S 



A /" 

d'^x£^^^= ^ (n + 1)(;S„,„ -|-am_i,„+i)f,„,„ (a_i,i = 0). (68) 

0^m+n^4 

Thus we have the relations between the coefficients of the Lagrangian and the equations of motion: 

am,n = (m + l)(a„i,ri + f3m+l,n-l) , (69) 

bm,n = {n+ l){l3rn.n + arn-l.n+l) ■ (70) 

As expected there exists an integrability condition relating the coefficients in the equation of motion. One 
can easily check that 

nam-i,n - mbm,n-i = . (71) 

The origin of this condition lies in the fact that TrSm-i.n — ^£m.n-i is a total derivative for n^m > 1. 
This means we have included a number of total derivatives in the Lagrangian that do not contribute to the 
dynamics. Indeed, at order m + n + I > 2 the total number of free parameters is not 2{m + n + 1) as one 
might naively expect, but m + n + 2 owing to the constraints given by Eq. 1711 

We now have the most general bi-galileon theory, given by Eqs. IHSl [51] and [551 Referring back to the 
discussion in the previous section, we note that only n couples directly to matter, through an interaction 
J d^x ttT. It follows that the equations of motion are given by 






0^m+n^4 
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where the am,n and the bm,n are related according to Eq. [7TJ 

We conclude this section by stating, without proof, the generalisation to N galileons. The analogue of 

X [(a^.a^^^TTi) . . . (a^„^9''"i^i)] . . . [{dp.d^^nN) . . . {dp^j^'^-TTN)] . (74) 

Therefore, the general iV-galileon Lagrangian is given by 



where 



r 



7 , L'mi,...,7nN 1 



?Tii + ...+?njY ^4 



•'-mi,...,mjv ~ (,'^mi,...,mjv^l + ■ • ■ + O^nii ,...,171^'^^ )^'>ni 



(75) 



(76) 



V. MAXIMALLY SYMMETRIC VACUA AND THEIR STABILITY 

Although a detailed study of bi-galileon phenomenology will be left to our companion paper [4l|, we 
will now initiate a study of maximally symmetric vacua. We will be particularly interested in establishing 
stability conditions associated with ghost-like excitations of the galileon fields about these vacua. 

We will allow for a non-zero vacuum energy so that the background stress energy tensor is given by 
r^ = —(7(5^. As discussed in section Hill a maximally symmetric solution has a background it field. 



7r(a;) 



-K'jyX ^X , O^Oi/TT 



''^ir'n^iv : 



(77) 



where k-^ = H^ — H^j^ is the difference between the Hubble parameters calculated in the modified gravity 
theory and in GR. For simplicity we assume that the background ^ field takes the same form 



t,(x) — k^Xf^x 



dudui 



' 2^^f^'' ' 



(78) 



for some constant k^. If our theory happens to be derived from a brane model, we expect this choice to 
correspond to choosing maximal symmetry in the bulk and on the branes. We now plug this into our field 
equations Eqs.[72]and[73l This yields 



E^{k^,k{) 
E^{k^,k{) 



-4cr 



E 



m^n 



4! 



(A — m — n)\ 



h.mi,n _ r. 



0^m+n^4 



m+n 



4! 



(4 — ?Ti — n)! 



h h™-k"- — n 



(79) 
(80) 



where we recall that «„.„ and b^^n are given by Eqs. l69l and l70l Generically these represent two independent 
equations in two unknowns. Being quartic polynomials, their solutions will typically contain sixteen distinct 
points in the real (A:,r,fc{) plane. Of course, there may be degenerate configurations with a continuum of 
solutions, or other configurations with no real solutions whatsoever. We will study these aspects in more 
detail in [4l[. For the moment, let us assume that a real solution exists. 

Can we really trust our solution? We can think of two things in particular that could spoil the background: 
gravitational back reaction and/or dynamical instabilities. The former refers to backreaction of the scalars 
onto the geometry whereas the latter refers to ghost-like excitations. To succinctly express suitable conditions 
on the background, it is convenient to consider what happens when we plug our background fields directly 
into the action. It turns out that 



(Tx (C. 



ff,C 



Ttt) 



- f / d^x Xfj.x^' j L{k,r, k^) , 



(81) 
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where we define the action polynomial, 

L{K,k^)^AaK- V K,,„fc,+/3™,„fc5) -- ■■ -rKk^. (82) 

^-^ V 2/ (4 — m — n ! ^ 

Now, it is easy to check, using Egs. 1691 and [701 that 

Or n 7" 

^^(/c^,fc5) = -— , E^{K,k^)^-—. (83) 

Of course this comes as no surprise since the equations of motion are designed to extremise the action. For 
the maximal symmetry case in question this simply corresponds to extremising L. 

Now let us return to the issue of backreaction of the scalars on to the geometry. This was discussed at 
the end of section IIIII Schematically, we expect 

TZyJm^^^]-^' E 0(l)(«™,n-3M2j^J°)C+'fc^" + 0(l)/3,„,„fc™fcr' (^4) 

0^m+n^4 

Recall that to justify neglecting the backreaction, we require Tg'^aiar ^ M'^i£h^„ ~ M^iH^. 

Assuming that we can indeed neglect the backreaction as described, we then need to ask if our background 
solution is stable against excitation of ghosts? To study the stability of a solution we need to examine its 
fluctuations, so we shift the fields so that they represent fluctuations on the appropriate background 

TT^Tf + TT, ^^C + C- (85) 

We claim that the background solution has a ghost-like instability if the Hessian of L{kT^, k^) 

/ d^L d^L \ 

Hess{L) = JL '""S't (86) 

\ dkfdk^ dkf / 

has a negative eigenvalue on the solution. Let us now prove this claim. 

The field re-definition (1851) takes us from the trivial to a non-trivial background, and gives rise to the 
following shift in the Lagrangian 

£^,^ -^ C^^^ + C'^^^ . (87) 

It is clear that the perturbed Lagrangian C'^ c has the same symmetries as /CTr,^, so we may once again write 

where 

'~'ni,n ~ \'^m,n'^ + Pm.n^J^m.n ■ (,0") 

Note that, due to the background equations of motion, linear terms in C^ ^ vanish; that is there is no 
contribution from m = n = 0. The equations of motion are given by 

E ^'■m,n£m,n ^ -ST , (90) 

/ , '^m,n^m,n = U j (,9iJ 



with the relations 



aL,n = {m + 1)«,„ + /3,'„+i,„-i) (/3,:,_i = 0) , (92) 

Kn,n = {n + l)(/3:„^„ + «™-i,„+i) (a'_i^, = 0) , (93) 
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and 



'^am-l,™ - "^Cn-l =0. 



(94) 



There is a linear map between the parameters in our original theory am,,m bm,n and the perturbed theory 
^'m m b'm n- ^Y direct substitution we find that for 1 ^ ?7i + n ^ 4, we have 



4 4 



4 4 



'^7n,n ~ / , / , "'^m,n 0,i^j , bm,n ~ / , / , ^^^m,n ' t>ij , 



(95) 



where 



^^'"■" - I 2) [mj [nj (4-«-j)! " '^« 



l^i-ni^J-" _ 



(96) 



Here the combination 



i\/m\{i — m)l. For A^ < 0, we extend the definition of "factorial" using 



the Gamma function, TV! = T{N + 1), and recall that l/r(iV + 1) = when iV is a negative integer. It 
immediately follows that we only get contributions from i + j ^ 4. 

To study the stability of the background, we consider the leading order piece of the perturbed Lagrangian 






total derivative - 



The background has a ghost-like excitations if 



K = 



2a\ 



10 



/3i 



10 



P'l 



10 



^01 



2/3^ 



01 



(97) 



(98) 



has a negative eigenvalue. However, using Egs. I92[ 1931 and the linear map Eg. I95[ we see that 



K 



'^10 '^01 



\ E 



1 



i+i 



4! 



2 -^ V 2/ (A-i-j)\ 



lOijk^ kc jbijk^kc 



(99) 



Differentiating Egs. [751 and [501 then using Eq. [531 one is able to show by direct comparison that 

1 



K 



dkTT dk^ 

dE(_ dE( 
dkTY dk^ 



-Hess{L) 



(100) 



Recall that the background solution contains a ghost if K has a negative eigenvalue, which is equivalent to 
Hess{L) also having a negative eigenvalue. Thus we have proven our claim*. 

We can understand this result geometrically. The function L{k.,r,k^) is a smooth polynomial function 
of two variables, whose extrema correspond to maximally symmetric vacua in our theory. Unstable vacua 
happen to be those for which the Hessian of L has a negative eigenvalue on the solution. This just means that 
the corresponding extrema are saddles and maxima. Or put another way, stable vacua are simply minima 
of the action polynomial, L. Note that a vanishing eigenvalue will generically lead to strong coupling at all 
scales as the corresponding scalar degrees of freedom will have a vanishing kinetic term. 

This geometrical insight enables us to make some generic statements about the spectrum of vacua in each 
case, using Morse theory. For example, if our galileon theory contains more than one vacuum solution, they 



The analogous result holds for N galileons. 
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cannot all be stable. This is because Morse theory tells us that you cannot have all extrema being minima 
if there is more than one extremum. This logic explains why the self-accelerating branch of the DGP model 
must have a ghost-like instability [T^ if the normal branch is stable. 

The action polynomial is an extremely useful construct. One can use it to build interesting theories 
from the bottom up, ensuring stability conditions from the outset. This will be the approach we take 
in our detailed study of phenomenology [4l|. Given a suitable action polynomial, we can reconstruct the 
corresponding action for the relevant theory by computing the following coefficients: 



4!m!n! dk^^ dk^ 



fe„=fce=0 



A:„=fef=0 



4a<5^"<5^' (101) 
(102) 



where m, n = 0, 1, 2, . . ., and we recall that we define a_i,„ — Pm,-! ~ 0- Note that since T:£m-i.n ^ S,£m,n-i 
is a total derivative for n,m > 1, we are free to set, say, /3i_„ = P2,n = P3,n = . . . = 0, without loss of 
generality. However, it is clear from Eq. llOll that this choice subsequently fixes am,n uniquely. 

The action polynomial can also be used to quickly calculate the equations of motion for fluctuations on 
a solution, given by Eqs. |90] and 1911 Given the action polynomial, we can compute the coefficients in the 
perturbation equations directly, using the following formulae, valid for 1 ^ m -I- n ^ 4 

"■" " ^ ' Almlnl dk^+'dk^ ^ ' 

b' - ( 2)"+" ('-"^-")' ^"'^"^^^ (104) 

'"■" " ^ ' 4!m!n! dk^dk^+' ^ ' 

As it is much easier to work at the level of the action polynomial, these formulae will prove invaluable to 
anyone wishing to study the phenomenology of maximally symmetric solutions and their fluctuations. 

VI. DISCUSSION 

We are now at the half way stage in our study of multiple galileon theory. We will continue with a detailed 
analysis of phenomenology in our companion paper |4l| . but for now let us take stock of where we are. 
We have generalised the original galileon theory [24| to include any number of coupled galileon fields. We 
have focused our attention on the case of two galileon fields, one of which is coupled to the trace of the 
energy-momentum tensor. We call this version of the theory, bi-galileon theory. Although the two galileons 
are coupled, perhaps even at quadratic order, we neglect any direct coupling between them and the graviton. 
This amounts to neglecting the backreaction of the scalars to background geometry. 

Through the coupling to the energy-momentum tensor, the galileon fields lead to a long range modification 
of the gravitational force. We have argued that co-dimension two braneworld models with an infra-red 
modification of gravity will typically admit a decoupling limit corresponding to a bi-galileon theory. Indeed, 
we have shown explicitly how the cascading cosmology model [3ll does exactly that. Of course, strictly 
speaking the cascading cosmology model is a 51? model with co-dimension one branes ISll . although it is 
expected to capture many of the salient features of the full 6D cascading gravity model [2^, |30| which does 
contain co-dimension two branes. Our interest in galileon modifications of gravity need not rely on hi ghe r 
dimensional models of gravity. We can, in principle, reconstruct our bi-galileon theory along the lines of [27| . 
and develop a fully covariant model of two scalar fields coupled to gravity and to matter in four dimensions. 
While the Galilean invariance will be broken by this procedure, we expect it to be restored in the limit 
Mpi — >■ oo, so we can still learn plenty about the model simply by looking at its bi-galileon limit. 

Our construction of the most general bi-galileon theory is very similar to the single galileon case, with 
a straightforward extension to N galileons. As was the case for the single galileon [24|, there are only a 
finite number of terms you can write down in the equations of motion. This is because requiring Galilean 
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symmetry, as well as no higher than second derivatives, forces us to build our equations of motion from 
scalar products of dfj^d^?: and dp^du^. These can be thought of as matrices from which we can only build a 
finite number of Cay ley invariants. 

As a first step towards a study of the phenomenology, we have considered maximally symmetric vacua 
and their stability. We have found a very elegant way of describing the physics in terms of a particular 
function, the action polynomial, which is closely related to the on-shell background Lagrangian. Stable 
vacua correspond to the minima of this polynomial. Computing its Hessian on a particular solution offers us 



a quick and easy stability test which we will make use of in the future [41| . We can also apply Morse theory 
to understand the generic spread of stable and unstable vacua. For example, it becomes almost trivial to 
see that if there is more than one vacuum, they cannot all be stable. 

Having motivated and constructed our general theory we conclude by referring the reader to our 
forthcoming companion paper [4l| for a detailed discussion of phenomenology. 



Note added: Whilst this paper was in its very final stages of preparation, 42| appeared, in which 
the galileon model is generalised to mixed combinations of arbitrary p- forms. It is clear that the results 
of section HVl of this paper can be derived from the formulae given in 42], for the case of multiple scalars. 
Both this work and our own might be considered to be the natural generalisations of |43|. 
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